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Abstract. After a short review of the Physics motivations for the study of open heavy flavour production 
in proton-proton, proton-nucleus and nucleus-nucleus collisions at the LHC, we present results on the 
expected performance of the ALICE experiment for charm and beauty production measurements. 
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1 Introduction 

The ALICE experiment will study nucleus-nucleus (AA) 
collisions at the LHC, with a centre-of-mass energy -^Snn — 
5.5 TeV per nucleon-nucleon (NN) collision for the Pb-Pb 
system, in order to investigate the properties of QCD mat- 
ter at energy densities of up to several hundred times the 
density of atomic nuclei. Under these conditions a decon- 
fined state of quarks and gluons is expected to be formed. 

The measurement of open charm and open beauty pro- 
duction allows to investigate the mechanisms of heavy- 
quark production, propagation and hadronization in the 
hot and dense medium formed in high-energy nucleus- 
nucleus collisions. The open charm and open beauty cross 
sections are also needed as a reference to measure the ef- 
fect of the transition to a deconfined phase on the pro- 
duction of quarkonia. Heavy-quark production measure- 
ments in proton-proton and proton-nucleus collisions at 
the LHC, besides providing the necessary baseline for the 
study of medium effects in nucleus-nucleus collisions, are 
interesting per se, as a test of QCD in a new energy do- 
main. 



2 Heavy-flavour production from pp to AA 



Heavy-quark pairs (QQ) are expected to be produced in 
primary partonic scatterings with large virtuality > 
(2mQ)^ and, thus, on small temporal and spatial scales, 
Z\t ~ Z\r ~ 1/g < 0.1 fm (for rUc = 1.2 GeV). In nucleus- 
nucleus reactions, this implies that the initial production 
process is not affected by the presence of the dense medium 
formed in the collision. Given the large virtualities, the 
baseline production cross sections in NN collisions can be 
calculated in the framework of coUinear factorization and 
perturbative QCD (pQCD) |2]. For the estimate of base- 
line production yields in nuclear collisions (to be used 
for performance studies and preparation of the analysis 



strategies), a scaling of the yields with the average num- 
ber (iVcou) of inelastic NN collisions (binary scaling) is 
usually assumed: 

d'^AA(pA)/dPtdj/ = (TVcoii) X d^iVpQ /dptdy . (1) 

The expected cc and bb production yields for differ- 
ent collision systems at the LHC are reported in the first 
line of Tabled 0. These numbers, assumed as the ALICE 
baseline, are obtained from pQCD calculations at NLO 0, 
including the nuclear modification of the parton distribu- 
tion functions (PDFs) (4j in the Pb nucleus (details on 
the choice of pQCD parameter values and PDF sets can 
be found in |3]). Note that the predicted yields have large 
uncertainties, of about a factor 2, estimated by varying 
the values of the calculation parameters. An illustration 
of the theoretical uncertainty band for the D meson cross 
section as a function of pt will be shown in section^ along 
with the expected sensitivity of the ALICE experiment. 

Several effects can determine the breakdown of binary 
scaling in pA and A A collisions. They are usually divided 
in two classes, that we discuss in the following. 

Initial-state effects, such as nuclear shadowing, the mod- 
ification of the parton distribution functions in the nucleus 
due to gluon recombination at small momentum fraction 
X. Initial-state effects can, at least in principle, be stud- 
ied by comparing proton-proton and proton-nucleus col- 
lisions. It has recently been argued that, indeed, at LHC 



Table 1. Expected QQ yields per event at the LHC, from NLO 
pQCD calculations ^. For p-Pb and Pb-Pb, PDF nuclear 
modiflcation is included and A^'coii scaling is assumed. 



colliding system 


PP 


p-Pb 


Pb-Pb 


-\/ Snn 


14 TeV 


8.8 TeV 


5.5 TeV 


centrality 




min. bias 


0-5% a'""' 


cc pairs 


0.16 


0.78 


115 


bb pairs 


0.0072 


0.029 


4.6 
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energy, gluon recombination may occur even in pp colli- 
sions and affect the charm production cross section j^. 

Final-state effects in AA collisions, due to the interac- 
tion of the produced partons with the medium. Partonic 
energy loss in the medium is the main example of such 
an effect. Believed to be at the origin of the jet quench- 
ing phenomena observed in Au-Au collisions at RHIC , 
energy loss is expected to depend on the properties of the 
medium (gluon density and volume) and on the proper- 
ties of the 'probe' parton (colour charge and mass). Due to 
the large values of their masses, charm and beauty quarks 
are qualitatively different probes with respect to light par- 
tons, since, on QCD grounds, the in-medium energy loss 
of massive partons is expected to be reduced relative to 
that of 'massless' partons (light quarks and gluons) 7,8 
IH]. In addition to that, since at LHC energy most of the 
measured light-flavour hadrons will originate from a gluon 
parent, heavy-flavour particles, such as D mesons, will pro- 
vide a tool to tag a quark parent. As pointed out in 10 , 
the comparison of the high-pt suppression for D mesons 
and for light-flavour hadrons should test the colour-charge 
dependence (quark parent vs. gluon parent) of parton en- 
ergy loss, while the comparison for B mesons and for light- 
flavour hadrons should test its mass dependence (massive 
parent vs. massless parent) — in section 0] and [S] we show 
some predictions from and compare them to the ex- 
pected ALICE sensitivity for these quenching studies. 



3 Heavy-flavour detection in ALICE 

The ALICE experimental setup [T1lll| was designed in 
order to allow the detection of D and B mesons in the 
high-multiplicity environment of central Pb-Pb collisions 
at LHC energy, where up to several thousand charged par- 
ticles might be produced per unit of rapidity. The heavy- 
flavour capability of the ALICE detector is provided by: 

— Tracking system; the Inner Tracking System (ITS), the 
Time Projection Chamber (TPC) and the Transition 
Radiation Detector (TRD), embedded in a magnetic 
field of 0.5 T, allow track reconstruction in the pseu- 
dorapidity range —0.9 < rj < 0.9 with a momentum 
resolution better than 2% for pt < 20 GeV/c and a 
transverse impact parameter^ resolution better than 
60 /im for pt > 1 GeV/c (the two innermost layers of 
the ITS are equipped with silicon pixel detectors)^. 

— Particle identification system; charged hadrons are sep- 
arated via dE/dx in the TPC and in the ITS and 
via time-of-flight measurement in the Time Of Flight 
(TOF) detector; electrons are separated from charged 
pions in the dedicated Transition Radiation Detector 

^ The transverse impact parameter, do, is defined as the dis- 
tance of closest approach of the track to the interaction vertex, 
in the plane transverse to the beam direction. 

^ Note that, for pp collisions, the impact parameter resolu- 
tion maybe slightly worse, due to the larger transverse size of 
the beam at the ALICE interaction point. This is taken into 
account in the studies presented in the following. 



(TRD) and in the TPC; muons are identified in 
the forward muon spectrometer covering in acceptance 
the range —4 < ry < —2.5. 

Detailed analyses ^31 , based on a full simulation of the 
detector and of the background sources, have shown that 
ALICE has a good potential to carry out a rich heavy- 
flavour Physics programme. In section^ we describe the 
expected performance for the exclusive reconstruction of 
D'^ — > K~7r+ decays in pp, p-Pb and Pb-Pb collisions, and 
the estimated sensitivity for the study of charm energy loss 
in Pb-Pb collisions. In sections El and we present the 
perspectives for the measurement of beauty production in 
central Pb-Pb collisions in the semi-electronic and semi- 
muonic decay channels. 

For all studies a multiphcity of dNch/dy = 6000 was 
assumed for central Pb-Pb collisions^. We report the re- 
sults corresponding to the expected statistics collected by 
ALICE per LHC year: 10^ central (0-5% cr'"'=') Pb-Pb events 
at rpb-Pb = 10^^ cm-^s-i and 10^ pp events at £^^^'^^ = 
5 X lO'^*' cm^^s^^, in the barrel detectors; the forward 
muon arm will collect about 40 times larger samples (i.e. 
4 X 10* central Pb-Pb events). 



4 Measurement of charm production and 
in-medium quenching 

One of the most promising channels for open charm detec- 
tion is the D" K^7r+ decay (and its charge conjugate) 
which has a branching ratio (BR) of about 3.8%. The ex- 
pected production yields (BR x dN/dy at y = 0) for D° 
(and D^) mesons decaying in a K^tt^ pair in central Pb- 
Pb (0-5% cr"''^ at = 5-5 TeV, in minimum-bias 
p-Pb collisions at ^Snn — 8.8 TeV and in pp collisions at 

= 14 TeV are, in the order, 5.3 x 10-\ 3.7 x 10"^ ^^d 
7.5 X 10^^ per event. 

Figure ^ (left) shows a sketch of the decay: the main 
feature of this topology is the presence of two tracks with 

This value of the multiplicity can be taken as a conservative 
assumption, since extrapolations based on RHIC data predict 
d7Vch/dj/ = 2000-3000. 




1.78 1.8 1.82 1.84 1.86 1.88 1.9 1.92 1.94 1.96 

invariant iUlass [GeV] 



Fig. 1. Schematic representation of the D K~n decay 
(left). Ktt invariant-mass distribution corresponding to lO'^ 
central Pb-Pb events (right); the background-subtracted dis- 
tribution is shown in the insert. 
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6 18 20 
[GeV/c] 

Fig. 2. pt- differential cross section per NN collision for D" 
production, as expected to be measured with lO'^ central 
Pb-Pb events 10* minimum- bias p-Pb events, and 10^ pp 
minimum-bias events. Statistical (inner bars) and quadratic 
sum of statistical and pt-dependent systematic errors (outer 
bars) are shown. A normalization error of 9% for Pb-Pb, 9% 
for p-Pb and 5% for pp is not shown. 
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Fig. 3. Sensitivity on d^a^ /dptdy, in pp at 14 TeV, com- 
pared to pQCD predictions obtained with different sets of in- 
put parameters: iric [GeV], the factorization and renormaliza- 
tion scales, in units of mt,c, and the PDF set. The comparison 
is shown as a 'data/theory' plot. Error bars are defined as in 
Fig. HI 



the nuclear modification factor of D mesons as a func- 
tion of pt 



impact parameters do 100 fiui. The detection strat- 
egy 23] to cope with the large combinatorial background 
from the underlying event is based on: 

1. selection of displaced- vertex topologies, i.e. two tracks 
with large impact parameters and small pointing an- 
gle Op between the D° momentum and flight-line (see 
sketch in Fig. ^ ; 

2. identification of the K track in the TOF detector; 

3. invariant-mass analysis (see pt-integrated distribution 
in Pb-Pb after selections in Fig. ^| . 

This strategy was optimized separately for pp, p-Pb and 
Pb-Pb collisions, as a function of the D° transverse mo- 
mentum |15II13| . As shown in Fig. [3 the accessible pt range 
is 1-20 GeV/c for Pb-Pb and 0.5-20 GeV/c for pp and p- 
Pb, with a statistical error better than 15-20% and a sys- 
tematic error (acceptance and efficiency corrections, cen- 
trality selection for Pb-Pb) better than 20%. More details 
are given in Ref. |15lll3j . 

For the case of pp collisions, the experimental errors 
on the pt-differential cross section are expected to be sig- 
nificantly smaller than the current theoretical uncertainty 
band from NLO pQCD calculations (estimated by varying 
the values of the charm quark mass and of the factoriza- 
tion and renormalization scales). The resulting 'data/theory' 
plot in Fig. 13 shows that this will allow us to perform a 
sensitive test of the pQCD predictions for charm produc- 
tion at LHC energy. 

We studied ^H] the sensitivity for a comparison of the 
energy loss of charm quarks and of massless partons by 
considering: 



^AA(Pt) ^ 



1 



/dpt 



coll/ 



dN^p/dpt 



(2) 



which is used to characterize the medium-induced high- 
pt suppression — in central Au-Au collisions at RHIC, 
Raa is found to be ~ 0.2 for both 7r° and charged 
hadrons for pt > 4 GeV/c 0; 

the heavy-to-light ratio of the nuclear modification fac- 
tors of D mesons and of charged hadrons: 



i?D/,(pt) = i??A(Pt)/i?Ubt). 



(3) 



In Fig.0]we compare our estimated sensitivity on i?AA 
and Ru/h to theoretical calculation results jTH] that im- 
plement radiative parton energy loss with medium den- 
sity described by transport coefficient values in the range, 
q = 25-100 GeV^/fm, expected for central Pb-Pb col- 
lisions at ^Snn = 5.5 TeV on the basis of quenching 
measurements at RHIC. The experimental uncertainties, 
reported in Fig. 0] for the case q = 50 GeVVfm and 
nic = 1.2 GeV, are discussed in detail in Refs. [TCltTHj . The 
effect of nuclear shadowing, introduced via the EKS98 pa- 
rameterization ytj , is clearly visible in the i?AA without en- 
ergy loss for pt^7 GeV/c. Above this region, only parton 
energy loss is expected to affect the nuclear modification 
factor of D mesons. The small difference between the theo- 
retical i?AA predictions for rric = and 1.2 GeV indicates 
that the charm quark behaves similarly to a light quark, 
as far as energy loss is concerned. Therefore, the enhance- 
ment of the heavy-to-light ratio R^/h is a sensitive mea- 
surement, free of mass effects, to study the colour-charge 
dependence of parton energy. As shown by the error bars 
in the figure, Rri/h can be measured with good accuracy 
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Fig. 4. Nuclear modification factor for D" mesons (top) and 
heavy-to-light ratio of the nuclear modification factors for D" 
mesons and for charged hadrons (bottom). Predictions \W} 
with and without the effect of the charm mass are shown for 
the medium density range q = 25-100 GeVVfm. Errors corre- 
sponding to the case 'g = 50 GeV'^/fm and mc = 1.2 GeV are 
reported: bars = statistical, shaded area = systematic. 



(as it is a double ratio (AA/pp) / (AA/pp), some common 
systematic uncertainties cancel out). 



5 Measurement of beauty production in the 
semi-electronic decay channel 

The production of open beauty can be studied by detect- 
ing the semi-electronic decays of beauty hadrons, mostly 
B mesons. Such decays have a branching ratio of ~ 10% 
(plus 10% from cascade decays b ^ c ^ e, that only pop- 
ulate the low-pt region in the electron spectrum) . The ex- 
pected yields (BR x dN/ dy at y = 0) for b — > e -1- X plus 
b ^ c (-> e + X) X' in central Pb-Pb(0-5% a'''"^) at 
^Snn = 5.5 TeV and in in pp collisions at -^/s — 14 TeV 
are 1.8 x 10^^ and 2.8 x lO^"' per event, respectively. 

The main sources of background electrons are: (a) de- 
cays of D mesons; (b) neutral pion Dalitz decays 7r° 
7e^e~ and decays of light mesons (e.g. p and lo); (c) con- 
versions of photons in the beam pipe or in the inner detec- 
tor layers and (d) pious misidentified as electrons. Given 
that electrons from beauty have average impact parame- 
ter do — 500 fim and a hard momentum spectrum, it is 
possible to obtain a high-purity sample with a strategy 
that relies on: 
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Fig. 5. Minimum-pt-differential production cross section per 
NN collision for B mesons at y — 0, as expected to be measured 
from semi-electronic decays with 10^ central Pb-Pb events. 
Statistical errors (inner bars) and quadratic sum of statistical 
and pt-dependent systematic errors (outer bars) are shown. A 
normalization error of 9% is not shown. Suppression predic- 
tions |10| with and without the effect of the beauty mass are 
shown for the medium density range q = 25-100 GeV^/fm. 



1. electron identification with a combined dE/dx (TPC) 
and transition radiation selection, which is expected to 
reduce the pion contamination by a factor 10^; 

2. impact parameter cut to reject misidentified tt* and 

from Dalitz decays and 7 conversions (the latter 
have small impact parameter for pt ~ 1 GeV/c); 

3. pt cut to reject electrons from charm decays. 

As an example, with do > 200 /im and pt > 2 GeV/c, the 
expected statistics of electrons from b decays is 8 x 10* 
for 10^ central Pb-Pb events, allowing the measurement of 
electron-level pt-differential cross section in the range 2 < 
Pt < 18 GeV/c. The residual contamination of about 10%, 
accumulated in the low-pt region, of electrons from prompt 
charm decays, from misidentified charged pions and 7- 
conversion electrons can be evaluated and subtracted us- 
ing a Monte Carlo simulation tuned to reproduce the mea- 
sured cross sections for pions and D° mesons. A Monte- 
Carlo-based procedure can then be used to compute, from 
the electron-level cross section, the B-level cross section 
d(T^(pt > Pt"")/dy QHI' In Fig. El we show this cross sec- 
tion for central Pb-Pb collisions with the estimated sta- 
tistical and systematic uncertainties. The covered range is 
2 < pf"' < 30 GeV/c. 

The predicted suppression of the B meson p™'"-diffe- 
rential cross section due to b quark energy loss is also plot- 
ted in Fig. El The transport coefficient rang e 25-100 GeVVfm 
is considered and the two bands represent the results for 
TOb = and 4.8 GeV; the two bands are well separated 
up to p""" ~ 15 GeV/c. The quenching predictions are 
shown only for illustration, since the study of the B me- 
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son suppression will have to be performed by using as a 
reference the cross section measured in pp collisions. The 
sensitivity of this study is currently being investigated. 

6 Measurement of beauty production in the 
semi-muonic decay channel 

Beauty production can be measured also in the ALICE 
forward muon spectrometer, —4 < rj < —2.5, analyzing 
the single-muon pt distribution and the opposite-sign di- 
muons invariant mass distribution |17II18| . 

The main backgrounds to the 'beauty muon' signal 
are tt^, K='= and charm decays. The cut pt > 1.5 GeV/c 
is applied to all reconstructed muons in order to increase 
the signal-to-background ratio. For the opposite-sign di- 
muons, the residual combinatorial background is subtracted 
using the technique of event-mixing and the resulting dis- 
tribution is subdivided into two samples: the low- mass 
region, M^+^- < 5 GeV, dominated by di-muons orig- 
inating from a single b quark decay through b c(— > 
fj-^)fj.^ (BDsamo), and the high-mass region, 5 < M^+^- < 
20 GeV, dominated by bb — > p~ jX^ , with each muon 
coming from a different quark in the pair (BBdiff). Both 
samples have a background from cc — )■ /i^/i^ and a fit is 
performed to extract the charm- and beauty-component 
yields. The single-muon pt distribution has three compo- 
nents with different slopes: K and tt, charm, and beauty 
decays. The first component is subtracted on the basis 
of the identified hadron spectra measured in the central 
barrel. Then, a fit technique allows to extract a pt dis- 
tribution of muons from beauty decays. A Monte Carlo 
procedure, similar to that used for semi-electronic decays, 
allows to extract B-level cross sections for the data sets 
(low-mass p^ n~ , high-mass p'^n~, and pt-binned single- 
muon distribution), each set covering a different B-meson 
Pt > region. The results using only the single muons 
are shown in Fig. |^ Since only minimal cuts are applied, 
the reported statistical errors (inner bars) are very small 
and the high-pt reach is excellent. The main sources of 
systematic errors (outer bars) are: corrections for accep- 
tance and efficiency, subtraction of the background muons 
from charged pion and kaon decays, and fit procedure to 
separate the beauty and charm components. 

7 Conclusions 

Heavy quarks, abundantly produced at LHC energies, will 
allow to address several physics issues, in pp, pA and AA 
collisions. In particular, they provide tools to: probe, via 
parton energy loss and its predicted colour-charge and 
mass dependences, the dense medium formed in Pb-Pb col- 
lisions; probe, in pp collisions, the pQCD calculations pa- 
rameters space; probe the small- a; regime of the PDFs, 
where saturation effects are expected to be important. 

The excellent tracking, vertexing and particle identifi- 
cation performance of ALICE will allow to fully explore 
this rich phenomenology, as we have shown with some spe- 
cific examples on D and B meson measurements. 




B p;"'" (GeV/c) 

Fig. 6. Minimum-pt-difi^erential production cross section per 
NN colhsion for B mesons with —4 < y < —2.5 in central 
Pb-Pb coUisions, as expected to be measured from the single- 
muon data set. Statistical errors (inner bars) corresponding to 
4 X 10* events and pt-dependent systematic errors (outer bars) 
are shown. A normahzation error of 10% is not shown. 
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